FUNCTION initParameters
  {
	/*
	This function is used to extract the specific parameter values
	from the init_bounded_number_vector to the specific variables
	used in the code.
	
	Note that you must call this routine before runnning the 
	simulation model to generate fake data.
	*/
	
	ro = mfexp(theta(1));
	dvariable h = theta(2);
	switch(int(cntrl(2)))
	{
		case 1:
			//Beverton-Holt model
			alpha = (4.*h/(1.-h));
			break;
		case 2:
			//Ricker model
			alpha = pow((5.*h),1.25);
		break;
	}
	
	
	//TODO Alternative parameterization using MSY and FMSY as leading parameters
	
	
	m = mfexp(theta(3));
	log_avgrec = theta(4);
	//log_recinit = theta(5);  //DEPRECATED
	
	prho=theta(5);
	varphi=theta(6);
	log_q = q_prior;

	if(verbose)cout<<"**** Ok after initParameters ****"<<endl;
	
  }
